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Growth of interfaces during vapor deposition is analyzed on a discrete 
lattice. It leads to finding distribution of local heights, measurable for any 
lattice model. Invariance in the change of this distribution in time is used 
to determine the finite size effects in various models The analysis is applied 
to the discrete linear growth equation and Kardar-Parisi-Zhang (KPZ) equa- 
tion. A new model is devised that shows early convergence to the KPZ dy- 
namics. Various known conservative and non conservative models are tested 
on a one dimensional substrate by comparing the growth results with the 
exact KPZ and linear growth equation results. The comparison helps in es- 
tablishing the condition that helps in determining the presence of finite size 
effect for the given model. The new model is used in 2+1 dimensions to 
predict close to the true value of roughness constant for KPZ equation. 
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Growth on a lattice from vapor can be represented in primarily two ways. 
It can be modeled as a lattice model where the atomic interactions are re- 
placed by simple growth rules [1, 2] , then obtain a growth equation based on 
various symmetries of the problem under consideration [1]. Other way is to 
construct the growth terms from the given growth rules for a lattice model 
at the coerce-grained time and length scales [3]. The KPZ equation was in- 
troduced to include lateral growth in growth equation [4]. It has attracted 
a lot of attention in the field of growth. There are many lattice models 
and numerical solutions claiming to belong to the same universality as KPZ 
equation[5, 6, 7, 8, 9, 10, 11]. In 1+1 dimensions, exponents can be exactly 
obtained[l]. However, in higher dimensions exact values are not obtained. 
Various lattice models and numerical solutions predict a range of values due 
to the finite size corrections. In the following we develop a method to deter- 
mine the existence of finite size effect in a model. A model that converges to 
its representative universality can be identified and hence correct exponents 
can be determined. 

A linear equation representing interface motion normal to the surface can 
be obtained in the frame of reference moving with the interface velocity by 
considering inter-planer hopping of ad atoms on the interface with a bias for 
downward or in-plane hopping toward step edge [3]. It has the form 



where, z/ explicitly depends upon F, and 77 is the noise due to the randomness 
in the deposition flux. It has the correlation given by < ?y(x, t)r](x.', t') >= 
2D5(x — x!)S(t — t'). The angular brackets denote the ensemble average of 
the contents. Eq. (1) is known as Edward- Wilkinson (EW) equation [12]. 
The lowest ordered non-linear correction to EW equation was introduced by 
Kardar , Parisi , and Zhang [4]. The resulting equation, 



is known as KPZ equation. This is a non-conservative equation. 

The steady state growth is characterized by roughness exponent a and z 
, determining the evolution correlations in time. One can measure a from 
the height-height (h-h) correlations, 
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where, correlation length ~ t 1 ^. In the limit x — > 0, / — >• 1. Thus for 
a large £(£), the plot of t) vs. a; on the log scale must be a straight line for 
small x on any scalable surface. Hence any lattice model should comply with 
this requirement for large enough length and time scales. Absence of straight 
region over large enough length and time scales for a lattice model indicates 
that the corresponding surface is not scalable and hence such a model cannot 
exactly follow the growth equation that it is supposed to represent. We 
elucidate this point in the case of models believed to represent EW and 
KPZ equations. Time exponent (3, where z = a//3 can be obtained by 
measuring the width over a substrate of length L as, w 2 = jjJ2x(h(x,t) — 
h) 2 = L 2a g (^y)> It can be shown that [1] for small times w 2 ~ t 2/3 . 

We first analyze the G(x, t) for some of the models representing KPZ 
and EW type growth . Finite size effects enter due to both , the substrate 
size L and the cutoff length a. The rules in a lattice model are sensitive to 
both these lengths. This results in to deviation of the growth dynamics from 
that of the growth equation, that the model is representing. Signature of this 
deviation is obtained in the non linear behavior of G(x, t) in the limit x « £ 
where £ is correlation length, from a log-log plot of G(x, t) Vs. x . We have 
measured this deviation by fitting straight lines on log-log plot of G(x, t) Vs. 
x for every interval of Ax =10. From these straight line fits y=mx + cq, 
we obtain Co as a function of average x. In the absence of a curvature in 
the log-log plot of G(x ) t) Vs. x , one must obtain c to be independent of 
average x. 

The models are briefly described below. In most of the known modelsfl], 
finite size effects are present. We have chosen the model introduced in ref- 
erence [5] along with a new model. This new model, we believe, converges 
to the KPZ dynamics during early growth as will be seen from the results 
ahead. 

a) KK model [5]: In this model growth proceeds by selecting a site ran- 
domly (this is the first step in all the models described here.). A particle is 



3 



accommodated at the site if the absolute height difference between the se- 
lected site after deposition and for each of the nearest neighbors is less than 
or equal to a number N. 

b) SC model: We introduce another SOS model which provides limited 
tunability with respect to the spread in the distribution. This helps in iden- 
tifying exponent values close to the true values in 2+1 dimensions. The de- 
position rules for the model are as follows. In 1+1 dimensions the deposited 
atom is accommodated if both its neighbors have at least same height as the 
deposited one. Otherwise, largest of the step differences at the site ,Sd, is 
obtained and accommodation is allowed according to the probability factor 
e -s d /(2a )_ j} ere a can b e varied as a tunable parameter. In 2+1 dimensions 
the deposited atom is accommodated if three or more neighbors have at least 
same height as its own. For other depositions the accommodation is decided 
from the largest of the four steps around the site using above exponential 
probability factor. Details of this model are described elsewhere [13] 

c) NN1 model [3] : This is a conservative SOS model. A particle after 
deposition is allowed to relax by hopping to a nearest neighbor site if it can 
lower its height. The hop is not allowed if the height of one or more of its 
nearest neighbors is equal or larger. 

d) HM model [14]: This conservative model is based on the models pro- 
posed in ref [14] . Here, in a growth equation that involves terms of the form 
V 2 f(x), the growth proceeds by allowing the particle to hop to the nearest 
site that has minimum value for f(x). Thus, f(x) is like a potential. For 
V 2 /i, f(x) = h(x). For, V 4 /i, f(x) = -V 2 h. 

Models (a),(b) are assumed to belong to KPZ universality and models 
(d) and (e) to EW universality. 

Fig. 1 shows the log-log plots of G(x, t) Vs. x for various models. The 
substrate lengths and the number of MLs is large enough to provide satu- 
ration of G(x,t) around x = 1000. c is measured between x =10 to 100. 
Fig. 2 shows plot of c for various models as a function of average x. As 
can be seen, KK and NN1 models do not show straight line behaviour in this 
length scales. SC and HM models are straight lines within the statistical 
error bars. Deviation from straight line behavior is the indication of non 
scalable dynamics of the growth due to finite size effects in these models. 
Clearly, exponents derived from models like HM or SC are reliable in the 
respective universality classes in 1+1 dimensions since , from the straight 
line behavior, these models follow the dynamics of the growth equation that 
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they represent. 

This method requires measurement of G(x,t) over large substrate length 
and large times. In higher dimensions it is increasingly difficult to perform 
such measurements. In order to facilitate the determination of finite size 
effects in higher dimensions, we introduce another measurement based upon 
the time invariance of distribution of height fluctuations. 

Consider a one dimensional scalable lattice with a lattice constant a. We 
define step at site i as 

5xi = hi - h i+ i (4) 

The local slope is then —5xi/a. Consider linear growth equation Eq. (1) 
in 1+1 dimensions. || — > (hi(t + At) — hi(t))/At on discretization. The 
r.h.s. is (8xi-i(t) — 5xi(t))/a + i]i(t). This relation predicts the value of h at 
(t + At) dependent on the discrete differences in steps at t. Let Ahi = hi(t + 
At) — hi(t). From the nature of the equation < Ahi >— 0. Further, function 
hi(t) = Hi(t) — H(t), where Hi(t) and H(t) are values of height and average 
height measured from the substrate. On the r.h.s. of the discrete linear 
equation, < 8xi-\ — 5xi >= and < rji >= 0. Thus the differential term and 
the noise term can be averaged to zero independently. In fact this is true 
for any conservative differential term. This observation is related to the fact 
that noise does not couple to conservative differential terms except with q = 
mode [15]. For KPZ equation r.h.s. is proportional to (5xi-i — 5xi + 5x 2 + rji). 
Between the terms {8xi-\— Sxi) and (Sx 2 ) , the latter term is b a times stronger 
where b is scaling parameter. Hence for large b, Ahi will be determined by 
(Sx 2 ). We will therefore consider only non linear term contributing to Ahi m 
KPZ equation. The difference between the conservative terms and the KPZ 
term is that, on the r.h.s. of the equation, < (5xi) 2 0. Since average over 
l.h.s. is zero, we must have < (Sxi) 2 +r]i >= 0. For a conservative term, the 
average over the term and the noise are independently zero. Hence, in the 
case of KPZ equation, the KPZ term couples with noise [2]. In order that 
the average on r.h.s. be zero, evaporation or vacancy addition is associated 
with the growth process . 

The distribution of Ahi i n the case of KPZ equation is determined by 
(Sxi) 2 + r]i. Since we are assuming a steady state growth, it is required that 
for the given time interval At, the distribution of Ahi, which is same as that 
of (Sxi) 2 + rji, must be independent of time. The time independence indi- 
cates that the random force as represented by the noise term is adequately 
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compensated for by the stabilizing growth term. The time dependence for 
the distribution will indicate that on the growing surface 1) the weightages 
of configurations are changing in time, and/or 2)new configurations are gen- 
erated as growth proceeds affecting the morphology on given scale. This 
will imply that the true steady state is not obtained. Normally one identifies 
steady state region by inspecting w vs. t on logarithmic scale. The beginning 
of linear region on this plot is considered as the onset of steady state region. 

In general, Ahi is a result of all the terms on the r.h.s. of a growth equa- 
tion. Hence, cross over regions will not be discriminated with respect to the 
distribution of Ahi. The distribution of Ahi at small and large length scales 
may differ. However, in either case it must be time independent irrespective 
of the dominant term at the given scale. 

Local configuration defining growth term is directly related to Ahi. This 
suggests that a measure of Ahi can be obtained by defining local height with 
respect to a local reference. Such height will respond to local changes in 
heights and help in providing a measure of Ahi. We define such a height as 
a height measured from average height of neighbors. (hi)i oca i = hi — (hi-i + 
h i+ i)/2, proportional to the difference between the local steps. Incidentally, 
the expression is similar to the second derivative of h(x) representing EW 
term. However, for any growth term, same expression can be used in defin- 
ing a measure of Ahi. The definition (hi)i oca i is more general. Note that 
{hi)iocai = (Sxi — 8xi-\)/2. In 1+1- dimensions, the (h-h) correlations in the 
discrete form are 

G(n) = ((hi - h l+n f) (5) 

We assume that the correlation length £ is very large compared to the length 
a. Using the definition of steps Eq.(4), 

G(2) =< 5x1 + 5x 2 i+l + 25x l 5x i+1 > . (6) 

Let < 8x1 >=< $ x 1+i >= $ 2 and , < 5xi5xi_i >= s5 2 where s is the coupling 
between the steps around a given site %. The distribution for 5xi is always 
symmetric around zero and time independent for an ensemble average. 

In the limit £ — > oo Eq. (3) reduces to G(x) = cx 2a where constant 
c — G(l). Hence Eq. (6) can be written as 

2 2Q = 2 + 2s (7) 

where G(l) = 5 2 in the discretized case. Coupling s uniquely determines a. 
Thus, for s=-l/2,0,and 1, a is 0, 0.5, and 1 respectively. This analysis can be 
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easily extended to higher dimensions. The relation between a and s remains 
unchanged over a square, cubic or hypercube lattice in higher dimension. 
We have independently verified it for the EW model. A rough surface will 
be characterized by some value of s. Expressing {hj)i oca i in terms of steps 
it can be shown that 4 < {hi)f ocal >= (2 — 2s)5 2 . This shows that the 
distribution of (hi)i oca i can be used as a measure of a and the definition of 
{hijiocai is applicable to any growth equation which has < a < 1. However, 
distribution of (hi)i oca i cannot be used to represent that of A/ij. Latter 
quantity is a result of local configuration dependent term and the fluctuation 
due to the noise term. Thus the appropriate measure of the distribution of 
Ahi will be, uncorrelated fluctuations in (hi)i oca i- Hence, we measure the 
distribution of change in {hi)i oca i = S(hi)i oca i over a time interval At > w(t) , 
ensuring the uncorrelated fluctuations. For any model expected to follow the 
KPZ , EW or any other conservative or non conservative growth equation, it 
is required that this distribution must be constant in time. 

In the following we apply this criterion of constancy to the same models 
described earlier. Although above discussion was for growth on one dimen- 
sional substrate , the corresponding criterion can be easily extended to higher 
dimensions. We have chosen At = 100 MLs in 1+1 dimensions and 60 MLs 
in 2+1 dimensions where At is the time difference considered in the sim- 
ulations . The time difference over which the constancy of distribution is 
tested is from 500 MLs to 5000 MLs in 1+1 as well as in 2+1 dimensions. 
In this time interval for all the models considered, ln(w) Vs ln{t) curve is 
linear implying a steady state growth region. In 2+1 dimensions, two sets 
of 6{hi) local are generated, one corresponding to x and other y direction, and 
added. The constancy of the distribution is checked by measuring the ra- 
tio of the values in the distribution at zero for 500 MLs and 5000 MLs i.e. 
P = iqq( J5oo-/5ooo ) ; where j t is the count of S(hi)i oca i at zero. The ratios P 

are obtained by averaging over large enough runs so that values of P are 
statistically discriminated. In 1+1 dimensions over 3000 runs are required for 
substrate size of L = 8000. We have also calculated sum of absolute values of 
p. = ioo( J '.5oo-^,5ooo ) measured between i=-4 to +4 values of the 5(hi) local as 

an additional measure of the constancy of the distribution of 5(hi)i oca i. Here, 
-fi,500 and -J^sooo are the counts at ith position in the distribution at 500 and 
5000 MLs respectively. This range (-4 to +4) is chosen because one of the 
models used in the present work provides values of 5{hi)i oca i only within this 
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range. This sum is denoted by P sum - 

In Table I we present the results from the measurement of distribution of 
$ (hi) local f° r different models in 1+1 dimensions. Along with P and P sum , we 
also measure a for each model. However, as is seen from Fig. 1, G(x,t) for 
KK and NN1 models are curved on the log-log scale. Hence, the measured 
a values are not reliable. Fig. 3 shows a plot of such a distribution for the 
model NN1 described above. 

KK Model: It shows a deviation of -0.5±0.06% and 6.86+1.0% in P and 
Psum respectively. For KK model the distribution becomes narrower in time 
(-ve Pq) indicating that number of tilted regions are growing with time. We 
have used N — 1 for KK model as height limitation [5]. 

SC Model: With a = 1.7 it shows significantly small spread compared to 
KK model in the distribution. We have observed that for other values of a, 
spread is larger. The advantage of this model is that it does not suffer from 
cross over effects. By varying a it is possible to get faster convergence to the 
KPZ dynamics. 

We have performed simulations between 200 MLs to 2000 MLs, 2000 MLs 
to 20000 MLs. Corresponding results have same trends i.e. for SC model the 
Po and Psum are less than 0.1% and 1.5% respectively while for KK model 
the absolute value is larger than 0.5% and 6%. Our analysis is based on the 
spatial discretization of the growth equation. It is known that [16] the real 
space group renormalization leads to different results than the KK result. 
Also for the KK model, in reference [17] the a value for 1+1- dimensions is 
0.489 although in [5] from saturation of width at larger L values it is given 
as 0.5. Our results show that the pathology associated with KK model is 
manifested in the form of change in the distribution of S(hi)i oca i in time. This 
in turn indicates that the finite size effects are dominant over the time scales 
considered. 

The results for conservative growth models also confirm to this behavior. 

HM Model and NN1 Model: HM model is like solving linear second or- 
dered growth equation locally. It is expected to follow the dynamics exactly. 
The Po and P sum values are indeed close to zero for this model. The NN1 
model restricts the minimization of h(x) due to the constraint that it is im- 
mobile if one or more nearest neighbors are present. This model has larger 
values of P and P sum . 

Above comparison between HM and NN1 models show that finite size 
effects do not allow convergence to EW dynamics in NN1 model due to the 
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constraint over the time scale studied. We consider the deviation from zero 
for P an d P SMm as the measure of convergence to the universality class that 
the model is belonging to. Larger the deviation farther is the model from 
the convergence. This convergence is with respect to the true dynamics of 
the growth equation, belonging to a particular universality class. Thus, SC 
model converged to the KPZ universality within the statistical error bar for 
the measurements of P and P sum , but KK model has not yet converged. 
Similarly, HM model has converged to EW universality but NN1 model has 
not. It is expected that asymptotically these convergences occur. However, 
within the time and length scales used for our measurements, scaling correc- 
tions do not allow true value measurements for KK and NN1 models in the 
KPZ and EW universality classes respectively. 

Above results assert that a model can be considered to have converged to 
its representative growth equation dynamics if the spread in its distribution 
of S(hi)i oca i is minimum. The advantage of this measurement is that it can be 
performed over a relatively small amount of growth compared to the G(x,t) 
measurement. Thus, even in higher dimensions the method can be used to 
determine the finite size effects in a model at an early stage. 

We have applied this method to determine an accurate value of a for the 
KPZ equation in 2+1 dimensions. For KK (N=l) model the P is 0.29+0.02% 
and P sum = 5.0 ± 1.0% with a = 0.402 ± 0.016. For SC model with a = 2.5, 
P is 0.02 ± 0.05% and P sum = 1.3 ± 1.1% with a = 0.355 ± 0.001. For KK 
model a is measured between x = 5 and 20 while for SC model it is between 
x = 2 and 50. We have plotted these results in Fig. 4. These results show 
that in 2+1 dimensions, true value of a for KPZ equation is close to 0.36. 

In conclusion, we have established a new criterion that can be applied to 
check the finite size effects at an early stage of growth. When applied to the 
models representing KPZ equation, it is seen that a new model , SC model 
in this connection can be adjusted to follow the KPZ dynamics accurately. 
This is true in both 1+1 and 2+1 dimensions. Based on this study it is 
seen that KK model predicts exponents that are away from the true values. 
Only those models that satisfy the condition of invariance of the distribution 
of Ah(x) will follow the dynamics of the representative growth equation 
correctly. These are the models that have already converged to the respective 
universality in terms of the underlying dynamics. We find that SC model 
converge to the universality much earlier. The KK model is not converged 
over the time and length scales used here. The HM model converges to EW 
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universality but NN1 does not. Hence, the exponent measurements based on 
converged models are reliable. 
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Figure 1: Plot of G(x,t) vs. x for KK model(N=l) (open circles), for HM 
model (open squares) , for NN1 model (filled circles), and for SC model 
(filled squares) in 1+1 dimensions. The curves are shifted along y-axis to 
avoid overlapping data. The growth is over 5.10 5 MLs with L=80000. 



Figure 2: Plot of Cq the straight line intercept on y- axis , as a function 
of average x. The plots are for the models in (1+1) dimensions, KK(N=1) 
(open squares), SC (+), NN1 (x) and, HM (filled squares). For the sake 
of comparison, values at x — 15 are adjusted to the same value for all the 
models. 



Figure 3: Plot of distribution of S(hi)i oca i for the model NN1 on semi log 
scale. The inner distribution (+) is for t = 500MLs and the outer one (x) is 
for t = 5000 MLs. 



Figure 4: Plot of G(x,t) vs. x for KK model(N=l) (open squares), and for 
SC model (filled squares) in 2+1 dimensions. 



Table 1: a values as obtained from height-height correlations, Pq the ratios 
of values , and sum of ratios P sum of 8{hi)i oca i at 500 MLs and 5000 MLs for 
different models in 1+1 dimensions. 
Model/Parameter a P in % P SU m m % 



KPZ and 



EW Equation 

KK(N=1) 

SC (a = 1.7) 

NN1 

HM 



0.5089 ± 0.012 -0.5 ± 0.06 6.8 ± 1.0 

0.5062 + 0.0015 0.07 + 0.02 1.1 + 0.4 

0.514 ± 0.02 1.82 ± 0.12 14.3 ± 2.0 

0.496 ± 0.002 0.06 ± 0.05 1.5 ± 1.2 
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